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ABSTRACT 


This  study  investigates  the  failure  modes,  failure  strengths,  and  failure 
criteria  of  laminated  composite  plates  with  stress  concentration  and  subjected 
to  bending  loading.  Graphite/epoxy  composites  are  used  for  the  present  study. 
Lamina  material  properties,  such  as  stiffness  and  strength,  of  the  composite 
are  obtained  by  experiments.  A  series  of  bending  tests  are  conducted  for 
laminated,  graphite/epoxy  composite  plates  with  and  without  a  hole  to 
investigate  their  failure  modes  and  strengths.  In  addition,  finite  element 
analyses  are  performed  to  compute  stress  distributions  around  holes  of  the 
composite  plates  subjected  to  bending  loading.  Based  on  the  stress 
computation,  a  couple  of  failure  criteria  are  examined  to  predict  the  failure 
strengths  of  composite  plates  with  stress  concentration. 
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B.  TRIANGULAR  PLATE  BENDING  ELEMENT 


L  INTRODUCTION 


Many  of  current  modem  technologies  require  materials  with  unusual 
combination  of  properties  that  cannot  be  met  by  conventional  metal  alloys, 
ceramics,  and  polymeric  materials.  This  is  especially  true  for  materials  that 
are  needed  for  aerospace,  underwater,  and  transportation  applications.  For 
example,  aircraft  engineers  are  increasingly  searching  for  structural  materials 
that  are  light,  strong,  stiff,  and  abrasion-,  impact-,  and  corrosion-  resistant. 
Rapid  increase  in  the  use  of  advanced  composites  is  primarily  due  to  their 
high  strength-to- weight  ratio  and  high  stiffness-to- weight  ratios,  making 
them  ideally  suited  for  the  high  technology  construction.  Naval  structure 
applications  of  composite  materials  are  typically  depending  on  this  weight 
saving  when  a  composite  material  is  used  to  fulfill  stiffness  or  strength 
requirements.  The  Navy  has  used  composite  materials  since  at  least  1946 
when  two  28  foot  personnel  boats  were  designed  using  a  laminated  plastic. 
Today,  composites  are  used  widely  in  the  fleet  and  include  items  such  as 
submarine  hatches  and  rudders,  MK  46  torpedo  propellers,  minesweeper 
hunters,  and  patrol  crafts  [Ref.  1]. 


Failure  of  engineering  materials  is  undesirable  for  several  reasons;  i.e, 
human  lives  that  are  put  in  jeopardy,  economic  losses,  and  the  interference 
with  the  availability  of  products  and  services.  Even  though  the  causes  of 
failure  and  behaviors  of  materials  may  be  known,  prevention  of  failures  is 
difficult  to  guarantee.  The  usual  causes  are  improper  material  selection  and 
processing,  inadequate  design  of  components,  and  their  misuse.  It  is  the 
responsibility  of  the  engineer  to  anticipate  and  plan  for  possible  failure  and, 
in  the  event  that  failure  does  occur,  to  assess  its  cause  and  then  take 
appropriate  preventative  measures  against  future  incidents. 

The  purpose  of  this  study  is  to  investigate  the  failure  modes,  failure 
strengths,  and  failure  criteria  of  laminated  composite  plates  with  stress 
concentration  and  subjected  to  bending  loading.  Both  experimental  and 
numerical  studies  are  conducted  for  the  study.  As  an  experimental  study, 
four-point  bending  tests  are  undertaken  to  obtain  lamina  material  properties 
as  well  as  to  examine  failure  modes  and  strengths  of  laminated  composite 
plates  with  or  without  holes.  As  a  numerical  study,  finite  element  analyses 
are  performed  to  determine  the  stress  distributions  around  the  holes.  A 
couple  of  failure  criteria  are  evaluated  using  both  the  numerical  and 
experimental  studies. 
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IL  LITERATURE  SURVEY 


In  recent  years,  there  have  been  significant  studies  to  determine  the 
stiffness  and  strength  of  laminated  composite  materials.  One  of  the  common 
experiments  used  to  define  the  responses  of  composite  materials  is  the 
flexure  test,  in  which  the  specimen  is  subjected  to  either  a  central  load  or 
two  quarter-point  loads.  In  many  cases,  results  of  these  tests  are  computed 
by  use  of  formulas  based  upon  a  homogeneous  material,  although 
bi-directional  (0°-90°)  specimens  are  employed,  hi  other  words,  the  geometry 
of  the  stacking  sequence  is  not  taken  into  account 

In  1949,  Hoff  [Ref.  2]  presented  flexure  formulas  which,  although  they 
were  based  on  an  elementary  strength  of  material  consideration,  appeared  to 
be  quite  adequate  to  describe  the  bending  of  unidirectional  or  bi-directional 
beams.  These  theoretical  equations,  in  conjunction  with  flexure  experiments, 
can  be  utilized  to  calculate  intrinsic  material  properties.  Certain  difficulties 
are  still  encountered  in  the  calculations  when  highly  anisotropic  composites 
are  tested.  Pagano  [Ref.  3]  reviewed  Hoffs  equations  and  expressed  them  in 
a  more  convenient  form  for  computational  purpose. 
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In  1974,  Haim  and  Tsai  [Ref.  4]  presented  graphs  to  determine  the 
laminate  stiffness  and  strength  for  the  [0°/+45°/90°]  boron  fiber/epoxy  and 
declaimed  that  a  set  of  generalized  graphs  for  all  practical  laminates  of  a 
given  composite  can  replace  the  current  limitation  of  one  set  of  graphs  for 
each  discrete  laminate. 

The  first  calculation  of  the  effect  of  a  circular  hole  on  the  stress 
concentration  in  a  plate  was  published  by  Kirsch  [Ref.  5]  in  1898.  He  dealt 
with  isotropic  plates  and  computed  the  stress  concentration  factor  of  3. 
Continuous  efforts  were  made  by  many  researchers  to  work  on  the  stress 
concentration  for  various  cases.  To  predict  the  strengths  of  laminated 
composites  containing  through  the  thickness  discontinuities,  two  related 
criteria  based  on  the  normal  stress  distribution  were  developed.  Whitney 
and  Nuismer  [Ref.  6]  considered  a  circular  hole  of  radius  R  in  an  isotropic 
plate  of  infinite  extent  with  the  origin  of  the  x-y  axes  system  at  the  center  of 
the  hole.  A  uniform  tensile  stress  was  applied  parallel  to  the  y-axis  at 
infinity.  As  the  first  failure  criterion,  they  assumed  that  failure  occurred 
when  the  stress  over  a  characteristic  distance,  dj,,  away  from  the 
discontinuity  was  equal  to  or  greater  than  the  strength  of  the  unnotched 
material.  They  proposed  the  second  failure  criterion  stating  that  failure 
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occurred  when  the  average  stress  over  a  distance,  a,,,  from  the  discontinuity 
equaled  to  the  unnotched  laminated  strength.  The  criteria  included  the  size 
effect  of  discontinuity.  Later  on,  Konish  and  Whitney  [Ref.  7]  found  that 
simply  scaling  the  stress  distribution  for  an  isotropic  plate  containing  a  hole 
by  a  ratio  of  the  isotropic  and  orthotropic  stress  concentration  factor  to  get  an 
approximate  stress  distribution  in  an  orthotiopic  plate  with  a  hole  might  lead 
to  significant  errors  even  close  to  the  hole  boundary,  particularly  if  the 
integral  of  the  stress  distribution  was  used  as  a  strength  criterion  in  the 
former  study.  They  proposed  a  simple  polynomial  which  was  obtained  by 
adding  sixth  and  eighth  order  terms  to  the  isotropic  solution  in  order  to  get 
an  approximation  for  the  normal  stress  distribution  adjacent  to  a  circular  hole 
in  an  infinite  orthotropic  plate. 

In  1981,  Hoff  and  Muser  calculated  the  stress  concentration  factors  for 
circular  holes  in  circular  composite  plates  of  cylindrical  orthotropy  when  the 
outer  edge  of  the  plates  were  subjected  to  a  uniform  uniaxial  traction.  They 
showed  that  the  stress  concentration  factor  can  be  plotted  in  a  diagram  as  a 
function  of  SB  (the  compliance  in  radial  direction)  and  S  (two  times 
tangential  compliance  plus  shear  compliance)  if  the  ratio  R  of  the  plate  to  the 
diameter  of  the  hole  was  prescribed.  They  also  made  a  conclusion  that  the 
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stress  concentration  factors  depended  primarily  on  the  local  arrangement  of 
the  fiber  at  the  ec^e  of  the  hole  [Ref.  8]. 

Chang  and  Chang  [Ref.  9]  presented  a  progressive  damage  model  to 
assess  damage  in  laminates  with  arbitrary  ply-orientations  and  to  predict 
ultimate  tensile  strengths  of  notched  laminates.  They  developed  a  nonlinear 
finite  element  program  for  laminates  containing  a  circular  hole.  They 
claimed  a  good  agreement  was  found  between  the  numerical  prediction  and 
the  experimental  data. 
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ffl.  EXPERIMENT 


This  section  explains  the  plate  bending  theory  for  laminated  composites 
to  obtain  composite  material  properties  from  bending  tests.  The  experimental 
apparatus  used,  the  laboratory  condition,  and  the  testing  procedure  are  also 
described. 

A.  PLATE  BENDING  THEORY 

Because  the  flexural  test  is  performed  to  obtain  material  properties  of  a 
graphite/epoxy  composite,  the  plate  bending  theory  for  laminated  composites 
is  reviewed  below. 

For  the  classical  plate  bending  theory,  the  following  assumptions  are 
made.  It  is  assumed  that  the  laminate  thickness  is  small  compared  to  its 
lateral  dimensions.  It  is  also  assumed  that  there  exists  a  perfect  bonding 
between  any  two  laminae.  That  being  so,  the  laminae  are  not  capable  of 
sliding  over  each  other  and  they  have  continuous  displacements  across  the 
interface.  It  is  made  another  important  assumption:  namely,  a  line  originally 
straight  and  perpendicular  to  the  laminate  midplane  remaines  so  after 
deformation. 
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Finally,  the  so-called  Kirchhoff  assumption  is  made  which  states  that 
in-plane  displacements  are  linear  functions  of  the  thickness  and  interlaminar 
shear  strains  in  the  x-z  and  y-z  planes  are  negligible.  These  assumptions 
reduce  the  laminate  behavior  to  a  two  dimensional  analysis  of  the  laminate 
midplane. 


The  three-dimensional  strain-displacement  relationships  are: 

_  _  du 
ex to 

C 

% 

dw 

E2‘& 

du  .  dv 


du  .dw 
„  dv  .dw 


.(i) 


Based  on  the  Kirchhoff s  assumption,  the  in-plane  displacements  are  linear 


functions  of  the  thickness  coordinate  z, 

“  “  “oM  +  ^1  (xO0; v  "  vo(x>>0 + . (2) 

where  Uq  and  v0  are  in-plane  displacements  of  the  midplane.  From  the 

Kirchhoffs  assumption,  Eq.  (3)  is  also  obtained. 

txz  m  F\ (x,y) +  zyz  "  F2(x,y)  +g£m°  . (3) 

It  follows  therefore  that  the  vertical  displacement  of  any  point  does  not 


change  in  the  thickness  direction.  Finally,  the  three  in-plane  strain 
components  are  expressed  as 


du  _  du0  d2w  0  .  . 

-  —  -  -r- -z — r-  -  Ev  +zkx 


dx  dx 


dv  ^0  d^w  0  .  . 


where 


OUq 

dx 

£vg 

dy 

top  f  dvp 
<)y  dx 


Equation  (4)  can  be  expressed  in  the  compact  form  as  given  bellow: 


ex  E*  kx 

Zy  “  ey  +2  ky 

e^J  l*’9' 


Now,  it  is  considered  a  composite  made  of  n  stacked  layers  or  plies.  Let 
t  be  the  total  thickness  of  the  laminated  composite  plate.  Then,  the 
constitutive  relationship  for  the  k-th  layer  is 

__  l°h "  [21*1®!*  . (6) 

where  (Q)fc  is  the  transformed  reduced  stiffness  matrix.  Substitution  of 
Eq.(5)  into  Eq.(6)  yields 

. (?) 

Since  a  laminated  composite  is  in  the  state  of  plane  stress,  there  are  three 
stress  components.  The  three  stress  components  produce  the  following  stress 


resultants: 
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(8) 


Nx  -f^oxdz 

Ny~f^nay<k 

N*y  m 

These  stress  resultants  have  the  dimensions  of  force  per  unit  length.  These 

resultants  give  the  total  force  per  unit  length  acting  at  the  midplane. 

Additionally,  moments  which  are  produced  by  the  stresses  with  respect  to 

the  midplane  are  also  applied  at  the  midplane,  and  they  are  written  as 
Mx  ~  °x2dz 

. 

Mxy  -  xxyxdz 

From  above  equations,  the  stress  resultants  are  written  as  a  summation  over 


n  plies,  i.e: 

Nx  ]  A\2  a16  c*  ^11  B\1  b\6  T  *x  ' 

Ny  -  A12  *22  A26  ey  +  B12  B22  B26  (10) 

Mxy  J  A16  A26  A66  e$y  B1 6  B26  B66  [kj^ 

or 

[M-Ulfe°] +  [*][*!  (11) 

Similarly, 

M  -  Plfe0]  +  tf>ll*l  (12) 

where 

Av’£l  (QiMk~zk-0 
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(13) 


Rewriting  Eqs.(ll)  and  (12)  together  gives 

[wl'fsolf  t  ] 

In  the  present  study,  four  different  laminates,  [0°]24,  r90°]24,  [0®/90®],2, 
and  [(+45°/-45°)]s,  subjected  to  pure  bending  moments  a  studied.  Since 
[0°]24  and  [90°]^  laminates  have  similar  stacking  properties,  they  will  be  put 
into  the  same  case  to  discuss: 

1.  Case  I  -  [0®]^  and  [90°]^  laminates: 

For  both  [0°]^  and  [90°]^  laminate  composite  subjected  to  pure 


bending  moments,  Eq.(12)  becomes: 

lAfl-[DPl  (15) 

because  the  matrix  [B]  will  become  zero  for  a  symmetric  lay-up. 

That  is, 

Mx  D\\  D\2  d\6 

My  -  Dyi  D 22  #26  (16) 

M*y  &26  d66  kxy 

Since  there  is  no  midplane  stretch, 

M— *W  (17) 


Substitution  of  Eq.  (17)  into  Eq.  (16)  and  inverse  of  the  resultant  expression 
gives 
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ex 

*>11  *>12  *>16 

-1 

Mx 

Zy 

-  -z 

*>12  D22  £>26 

My 

ex y 

D16  D26  D66 

Mxy 

Because  D16  and  D*  are  equal  to  zero  for  unidirectional  laminates  (0° 
or  90°)  and  cross-plied  laminates  (0°/90°),  Eq.  (18)  becomes: 


ex 

'*>11 

*>12  0 

-1 

'  Mx  ' 

Zy 

-  -z 

d12 

D2  2  0 

My 

exy 

0 

0  *>66 

Mxy 

where: 


A-det 


\Dh 

Dh 

0 

’  Mx  ' 

m  z 
A 

D12 

*>22 

0 

My 

0 

0 

*>66 

Mxy 

*>11 

*>12 

0 

*>12 

*>22 

0 

0 

0 

*>66 

D* 


u 


^22^66 


.(19) 


D*j2*  -D,2Da 


D*22*  ^11^66 


^11^22'^212 


•(20) 


That  is. 


ex 

ty 

zxy 


*>il*>22*>66 


'22"66 
12 
0 


"*>12*>66  *>11*>66 


Uu66 
L* 

0 


0 

0 

*>1  1*>22  “*>I2 


Mx 

My 

Mxy 


•(21) 
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If  subjected  to  a  uniaxial  pure  bending  moment:  Mx-M0;  the 

strains  are  expressed  as  functions  of  the  layer's  thickness,  material 


properties,  and  the  bending  moment,  hi  other  words, 

«*  ■  -TT-Wo  . (22) 

D  11 

For  the  four-point  bending  test  as  shown  in  Fig.  1, 

wo-f; 

where : 

p*  load  (lb) 

a*  length  of  moment  arm  (in) 
w-  plate  width. 


Rewriting  Eq.(22)  for  D„  gives 

D11 


zPa 

2e*w 


(23) 


For  [0°]^  and  [90%,  laminates,  Dn-2*h3*Qn/3  ,  D12«2*h3*Q12/3  ,  and 
D22*2*h3*Q2a^3  where  h  is  a  half  of  the  total  plate  thickness.  As  a  result: 


Q  ii- 


JSl _ 2 gz,£. 

Ah^zxMf  4 h?w  Zx 


(24) 


Since  oy-Txy-  0.  Eq.  (6)  yields 


3 pea 

ox  - 

Ahrw 


(25) 
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Figure  1. 4-point  bending  test 
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2.  Case  II  -  [0°/90°]u  cross-ply  laminate  : 

The  derivation  is  similar  to  that  for  [0°]  and  [90°]  laminates. 

However,  since  [0°/90°]12  is  nc  .  .symmetric,  the  matrix  [B]  will  not  reduce  to 
zero  and  the  situation  becomes  more  complicated 

The  stacking  sequence  and  the  thickness  of  each  corresponding  ply 
are  shown  in  Fig.  2.  From  Equation  (13),  the  [A],  [B],  and  [D]  matrices  for 
[0°/90°]12  are  obtained  as  when  each  layer  has  the  same  thickness. 
A„-h*[(Qu)  +  (Qja)],  A22*A11 

AI2«2h*(Q12),  Aw*A26*0 

A„-2h*(Q J, 

Bu-hJ*[(Q„)  +4*  (QJ1/96,  Bjj-h^CQ,,)  +  (QJ1/96 

B,2-5h2*(Q12)/96,  B1#-B*-0  (26) 

B«-5h2*(QJ/96, 

D„*h3*[(Q„)  +  (QnJl/3,  DB>Dn 

Da-2^QJ/3,  D„-Dm-0 

D«-2hJ*(QJ/3  . 

3.  Case  III-  [+45V-45*]  laminates : 
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For  the  [(+45°/-45°)6]j  laminated  composite  plate,  the  derivation 
procedure  is  also  more  complicated  than  the  [0°  ]  and  [90®]  laminates.  Due 
to  symmetry,  matrix  [B]  becomes  null  and  the  detail  of  matrix  [D]  is 
provided  below: 

D*n  0*12  0*16  lr  Mx  ' 

D*i2  0* 22  0*26  . (27> 

D*  16  ®*66  Mxy 

where : 


D*„“D22D66-D226, 

D*,2*-D12D66+D16D26, 

^*16m^12^26'^16^22» 

D*22=Dj  p^l6t 

D*26--D.iD*+D12D16, 

and 

Dn-h’m,)  +  (QJ  +2*(Q12)  +4*(QJ]/6, 

^22*^11* 

DI2-h3*[(Qn)  +  (Qjj)  +2*(Q12)  -4*(QJ]/6,  . (28) 

Di«“D26»h3*[(QII)  -  (QJ1/48, 

D««h3*[(Q„)  +  (QJ  -2*(Q12)  ]  /6 . 

In  this  thesis,  Eq.  (24)  and  the  basic  equations  of  material  properties  are  used 
to  to  compute  the  elastic  modulii  E„,  E^,  ....etc.  The  details  to  calculate 
these  elastic  modulii  are  given  in  the  later  chapter. 
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B.  APPARATUS 


All  the  bending  tests  were  conducted  in  a  laboratory  whose  ambient 
temperature  is  22.5+2.0°C  with  an  average  relative  humidity  equal  to 
41+5%.  The  MTS  810  Material  Test  System  was  used  to  conduct  all  the 
bending  tests. 

The  materials  used  during  the  experimental  testing  were  graphite 
fiber/epoxy  composites  with  [0°]^,  [90°]^,  [0°/90°],2,  and  [(+45°/-45°)6]s 
laminates.  The  experimental  samples  were  cut  from  different  orientation 
laminate  plates  with  nominal  dimensions  as  12"  by  12"  square  and  0.15" 
thick  into  small  plates  with  length  of  12"  and  width  of  1".  Some  samples 
were  drilled  a  quarter  inch  hole  at  the  center.  A  fixture  was  designed  to  hold 
the  composite  samples  in  a  simply  supported  configuration.  Each  sample 
was  positioned  on  the  top  of  the  supporters.  A  four-point  bending  test  was 
performed  as  shown  in  Fig  1. 

Several  strain  gauges  (1/4  inch  CEA-  13-250UN-350,  with  gauge  factor 
2.120  +  0.5%)  were  mounted  on  both  sides  of  the  samples.  One  strain  gauge 
was  placed  at  the  center  of  each  side  of  [0°],  [(0°/90°)Js,  and  [+45°/-45°)6]s 
laminae  (without  hole)  for  the  purpose  of  obtaining  the  properties  of 
laminated  composite.  A  [90°]  laminate  sample  (without  hole)  was  mounted 


with  two  strain  gauges  on  each  side,  one  in  the  longitudinal  direction,  and 
the  other  in  the  transverse  direction. 


Additional  strain  gauges  were  mounted  on  each  sample  with  a  1/4"  hole 
at  the  center  to  measure  the  strain  response  of  the  sample  under  load  more 
accurately.  Figures  3  to  6  show  the  locations  of  the  strain  gauges  on  each 
sample.  Tests  of  these  samples  were  conducted  with  the  MTS  material 
testing  machine  at  a  crosshead  speed  of  0.12  in/min  ( or  0.305  mm/min).  The 
MTS  machine  provided  readings  and  a  force-displacement  print-out  on  the 
x-y  plotter  (HP  7054A  x-y  Recorder)  for  the  applied  force  and  the 
measurement  of  the  deflection  of  the  sample.  Strain  gauge  outputs  were 
connected  to  a  Measurements  Group  SB- 10  Switch  &  Balance  Unit,  and 
readouts,  in  microstrain,  provided  by  Measurement  Group  P-3500  Strain 
Indicator.  The  results  were  read  and  recorded  manually  upto  the  fracture  of 
the  sample. 


L. 
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on  the  compression  side: 


on  the  tension  side: 


Figure  4.  Strain  gauge  locations  for  [90]  laminate 
without  a  hole 


Figure  5.  Strain  gauge  locations  for  [0]  and  [0/90] 
laminates  with  a  1/4"  hole  at  the  center 
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IV.  EXPERIMENTAL  RESULTS 


This  section  presents  the  results  as  obtained  from  the  experiments. 
While  some  description  of  the  experimental  and  calculated  results  is 
provided  here,  a  more  detailed  explanation  and  physical  interpretation  of  the 
results  will  be  given  in  the  later  chapter.  In  order  to  ensure  consistency  in 
recording  loading  responses,  samples  were  placed  as  uniformly  as  possible 
with  strain  gauges  in  the  same  positions.  The  same  crosshead  speed  and  the 
same  procedures  were  also  employed  for  each  separate  test.  The  results  are 
presented  graphically  and  in  a  tabular  form. 

A  SPECIMENS  WITFOUT  A  HOLE 

The  four  different  laminated  composite  samples  without  a  hole  were 
tested  to  provide  a  basis  to  compute  the  material  properties  and  to  compare 
with  samples  with  a  hole.  Tables  1  to  5  show  the  loads  and  strains  recorded 
from  each  experiment  for  the  samples.  The  load  vs.  strain  relationship  is 
plotted  in  Figures  7  to  14. 
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Table  1.  results  for  a  [0]  laminate  without  a  hole 


load  (lb) 

gauge#  1** 

gauge  #  2 

55.5 

395 

-399 

104 

736 

-745 

150 

1,015 

-1,057 

256 

1,760 

-1,795 

298 

2,035 

-2,080 

354 

2,396 

-2,456 

402.5 

2,712 

-2,786 

434 

2,916 

-3,001 

497 

3,319 

-3,419 

525.5 

3,506 

-3,617 

596.5 

3,950 

-4,089 

627 

4,138 

-4,289 

694 

4,561 

-4,742 

717 

4,701 

-4,892 

780 

5,082 

-5,303 

811.5 

5,271 

-5,509 

887 

5,729 

-6,009 

933 

5,998 

-6,305 

986 

6,315 

-6,654 

1,040 

6,650 

-7,007 

1,070 

6,836 

-7,193 

1,100 

6,995 

-7,366 

1,130 

7,177 

-7,565 

1,166 

7,385 

-7,779 

1,200 

7,558 

-7,864 

*:  failure  occurred  at  approximate  1210  lbs. 


**:  see  figure  3  for  the  strain  gauge  locations. 
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Table  2.  results  for  a  [90]  laminate  without  a  hole 


load(  lb ) 

gauge#  1** 

gauge  #2 

gauge  #3 

gauge  #4 

10 

22 

-941 

-18 

996 

20 

55 

-2,264 

-46 

2,438 

25 

70 

-2,924 

-62 

3,160 

* 

*:  failure  occurred  at  approximate  28  lbs. 


**:see  figure  4  for  the  strain  gauge  locations. 


Table  3.  results  for  [(0/90)]  ,2  laminate  without  a  hole 


ioad(Ib) 

gauge#  1** 

gauge  #2 

29 

-798 

755 

59 

-1,644 

1,540 

90 

-2,518 

2,344 

119.5 

-3,332 

3,078 

149 

-4,158 

3,835 

180 

-5,015 

4,601 

208 

-5,822 

5,317 

240 

-6,750 

6,126 

267 

-7,666 

6,914 

297 

-8,652 

7,749 

326 

-9,756 

8,660 

361 

-11,334 

9,944 

391 

-12,962 

11,255 

* 

*:  failure  occurred  at  approximate  401  lbs. 


**:  see  figure  3  for  the  strain  gauge  locations. 
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Table  4.  results  for  [+45/-45]-(l)  laminate  without  a  hole 


load  (lb) 

gauge#  1** 

gauge  #  2 

30 

1,767 

-1,803 

60 

4,678 

-4,871 

90 

8,654 

-9,388 

107 

13,759 

-16,064 

* 

*:  no  breakage  occurred,  quit  recording  due  to  a  very  large  deflection. 


**:  see .  gure  3  for  the  strain  gauge  locations. 


Table  5.  results  for  [+45/-45]-(2)  laminate  without  a  hole 


load  (lb) 

gauge#  1** 

gauge  #2 

4.5 

339 

0 

29.5 

1,544 

-1,272 

61 

3,099 

-2,969 

90 

4,576 

-4,633 

122 

6,425 

-6,805 

150 

8,057 

-8,799 

180 

9,767 

-10,980 

210 

11,795 

-13,708 

240 

14,783 

-18,065 

* 

*:  no  breakage  occurred,  quit  recording  due  to  a  very  large  deflection. 


**:  see  figure  3  for  the  strain  gauge  locations. 
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B.  SAMPLES  WITH  A  1/4-INCH  HOLE  AT  THE  CENTER 


Similarly,  three  different  laminated  composite  samples,  [0°],  [(0°/90°)]12, 
and  [(+45°/-45®)]I2  with  a  quarter  inch  hole  at  the  center  of  the  plate  were 
tested.  Five  strain  gauges  were  mounted  with  three  on  the  tension  side  and 
two  on  the  compression  side  for  each  sample.  The  three  gauges  on  the 
tension  side  were  placed  with  two  at  a  distance  of  O.S  in.  away  from  the 
center  and  another  one  at  1.25  in.  The  locations  of  strain  gauges  are  shown 
in  Figures  5  and  6.  The  results  for  each  sample  were  recorded  in  order  to 
compare  them  with  the  data  obtained  from  the  previous  section  (section  A) . 


Table  6.  results  for  a  [( +45/-4S)  12] s  laminate  with  a  1/4-inch 
hole  at  the  center 


load(Ib) 

gauge#  1** 

gauge  #2 

gauge  #3 

gauge  #4 

gauge#  5 

0 

0 

0 

172 

0 

29.5 

-2,910 

-2,834 

2,997 

2,683 

2,999 

49.5 

-5,151 

-5,022 

5,210 

4,711 

5,369 

69.5 

-7,834 

-7,686 

7,908 

7,066 

8,180 

84.5 

-10,783 

-10,670 

10,857 

9,465 

11,247 

90 

•12,282 

-12,151 

12,320 

10,619 

12,748 

* 

*:  no  breakage  occurred,  quit  recording  due  to  a  very  large  deflection. 


**:  see  figure  6  for  the  strain  gauge  locations. 
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Table  7.  results  for  a  [0]  laminate  with  a  1/4-inch  hole  at  the  cento* 


load  Ob) 

gauge#  1** 

gauge  #2 

gauge  #3 

gauge  #4 

gauge  #  5 

0 

0 

0 

167 

0 

0 

30 

-405 

-428 

578 

444 

425 

49.5 

-703 

-736 

868 

753 

720 

72 

-1,049 

-1,099 

1,210 

1,116 

1,068 

91 

-1,341 

-1,403 

1,496 

1,418 

1,358 

112 

-1,656 

-1,726 

mnm 

1,737 

1,666 

133 

-1,996 

-2,070 

2,127 

2,075 

1,989 

152 

-2,277 

-2,352 

2,398 

2,353 

2,254 

171 

-2,579 

-2,649 

2,688 

2,648 

2,532 

190.5 

-2,882 

-2,949 

2,977 

2,943 

2,813 

210.5 

-3,190 

-3,252 

3,270 

3,240 

3,095 

230 

-3,509 

-3,575 

3,542 

3,568 

3,366 

250 

-3,822 

-4,133 

3,838 

4,164 

3,574 

270 

-4,193 

-4,530 

4,187 

4,600 

3,827 

290 

-4,520 

-4,878 

4,480 

4,954 

4,057 

310 

-4,853 

-5,227 

4,775 

5,320 

*4173 

320 

-5,412 

-5,420 

5,551 

5,500 

*2655 

*:  gauge  #  5  indicates  failure  between  310  and  320  lbs. 


**:  see  figure  5  for  the  strain  gauge  locations. 


Table  8.  results  for  a  [0/90]  laminate  with  a  1/4- inch  hole 
at  the  center 


load  (lb) 

gauge#  1** 

gauge  #2 

gauge  #3 

gauge#  4 

gauge  #  5 

20.5 

-552 

-538 

527 

542 

496 

41 

-1,110 

-1,087 

1,044 

1,072 

979 

60 

-1,627 

-1,594 

1,519 

1,557 

1,418 

80.5 

-2,184 

-2,133 

2,026 

2,072 

1,882 

100.5 

-2,741 

-2,670 

2,533 

2,585 

2,339 

131 

-3,567 

-3,462 

3,286 

3,342 

3,014 

160.5 

-4,363 

-4,239 

4,075 

3,679 

190 

-5,165 

-5,038 

4,734 

4,826 

4,360 

220 

-5,998 

-5,843 

5,486 

5,563 

5,047 

* 

*:  failure  occurred  at  approximate  228  lbs. 


**:  see  figure  5  for  the  strain  gauge  locations. 


C.  CALCULATIONS  FOR  MATERIAL  PROPERTIES 
1 .  Poisson's  ratio  vM: 

By  definition,  Vg  is  the  Poisson's  ratio  for  the  transverse  strain  in  the 
j-direction  when  stressed  in  the  i-direction,  which  means 
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for  a  stress  applied  in  the  2-direction.  Applying  the  result  of  the  [90°] 
laminate  experiment  to  Equation(29),  the  average  value  of  v2I  is  found  to  be 
0.02225. 


2.  Elastic  modulii  En  and  E^,  &  Poisson's  ratio  v12: 

Here,  in  order  to  obtain  the  material  properties,  Eu  and  E^,  a 

regression  analysis  using  a  FORTRAN  program  was  performed  to  linearize 

the  ratios  of  load  vs.  strain  for  [0°]  and  [90°]  laminates.  As  a  result, 

Cp/cjc)Qo  "  ±154.95  x  103  and  0  "  ±7  245  x  103 

Since 

En 

Ql1  ’  l~vl2v2l 

q  -  EZ±  . (30) 

22  1‘v12v21 
n  v12E22 
QU  3—v12v21 

plugging  Eq.  (24)  and  Eq.  (30)  into  Eq.  (29)  results  in  the  following 
relationship: 

£‘l£222  -  ±(.-2s~Wu)0o  . (31) 

^22~^’llv21  4h*/> 

After  calculations,  results  for  these  two  elastic  modulii  are: 


Eu-  6.273  x  10*  psi  ;  E*-  0.3649  x  106  psi 


From 

E\\ 

vI2"£TTv21 

^22 

we  also  get 

v12-0.3825. 

3.  Shear  modulus  G12: 

A  Matlab  program  was  written  in  order  to  compute  the  material 
property,  G12.  The  experimental  result  of  load  vs.  strain  ratio  of  the 
[(+45°/-45°)Js  laminate  was  compared  with  the  load  vs.  strain  ratio 
calculated  from  the  program  using  a  guessed  value  of  G,2.  The  Incremental 
Search  method  was  employed  to  determine  the  correct  guess  value  for  GI2 
iteratively.  The  shear  modulus  G,2  obtained  from  this  way  was  1.045xl06psi. 


D.  CALCULATIONS  FOR  STRESS  &  STRESS  CONCENTRATION 
FACTOR  (SCF) 

The  failure  load  for  each  sample  is  shown  in  Table  9.  Here,  stresses  are 
calculated  at  z*  +h  or  -h  (at  the  top  or  bottom  of  the  sample)  for  each 
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sample.  From  Eq.  (25),  the  failure  stress  for  [0°]24,  [90°]24,  and  [(+45°/-45°)Js 


laminae  are 

a0/wo  ~  1  840x10 5  psi/per  unit  length, 
a0/w  -  1.26x10s  psi/per  unit  length, 
a 90/wo  “  5.08xl03  psi/per  unit  length, 
a45/wo  ~  61  86x103  psi/per  unit  length, 
a45fw  “  29-67*103  psi/per  unit  length. 


where  subscripts  wo  and  w  indicate  "without  a  hole"  and  "with  a  hole", 
respectively. 


Table  9.  dimensions  and  failure  loads  of  samples 


laminate 

length 

failure  (lb) 

[0]„„ 

5" 

1.5" 

1,210 

[0L 

10" 

3" 

310 

■H 

8" 

2" 

28 

■BH 

10" 

3" 

401 

■H 

10" 

3" 

228 

■UiliLIH 

10" 

3” 

*240 

■SSEEH 

10" 

3" 

*90 

*:  failure  due  to  a  very  large  deflection. 
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For  [0°/90°],2,  the  matrices  [A],  [B],  and  [D]  were  calculated  from  Eq.(26). 
Plugging  these  data  into  Eq.(14)  and  Eq.(7),  the  failure  stress  for  the 
[0°/90°]l2  laminate  are 

a0/90 -iop/wo  ■  213.86rl03  psi /  per  unit  length, 

<*0/90 -bottom/wo  “  12.53x103  psi/  per  unit  length, 
a  0/90- top/w  “  162.13x103  psi/per  unit  length, 

<*0/90 -bottomfw  “  9.49x103  psi/per  unit  length. 

Based  on  these  data,  the  experimental  stress  concentration  factors  (SC F)  for 
laminated  composites  are  1.46  for  the  unidirectional  plate  and  1.32  for  the 
cross-ply  plate  if  the  failure  occurs  due  to  the  largest  stress  at  an  edge  of  the 
hole. 
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load  lb 


Figure  7.  Stress-  strain  diagram  for  [0]  laminate  without  a  hole 

O  w 
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load- (lb) 


load  lb 


strain  um 


Figure  9.  Stress-  strain  diagram  for  [90]  laminate  without  a  hole 
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load  to 


Figure  1 1 .  Stress-  strain  diagram  for  [0/90]  laminate 
with  a  1/4"  hole  at  the  center 
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toad  tb 


-2  -1.5  -1  -0.5  0  0.5 

strain  um 


1.5 
x  10* 


Figure  13.  Stress-  strain  diagram  for  [+45/-45]  laminate  without  a  hole 
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strain  um 


Figure  14.  Stress-  strain  diagram  for  [+45/-45]  laminate 
with  a  1/4"  hole  at  the  center 


V.  FINITE  ELEMENT  ANALYSIS 


A.  OVERVIEW 

Since  the  earliest  development  of  the  finite  element  method,  a  large  of 
{date  bending  elements  have  been  developed  and  reported  in  the  literature. 
For  example,  Batoz  et  al.  [  Ref.  10],  with  the  purpose  of  identifying  the  most 
effective  element  for  the  thin  plate  analysis,  presented  an  assessment  of  flat 
triangular  plate  bending  elements  with  displacement  degrees-of-freedom  at 
the  three  corner  nodes  only.  They  claimed  that  the  most  efficient  and  reliable 
three-node  plate  bending  elements  were  based  on  the  discrete  Kirchhoff 
theory  and  a  hybrid  stress  model.  Later  on,  he  [Ref.  11]  also  presented  the 
explicit  expression  of  the  stiffness  matrix  of  the  nine  degree-of-freedom 
plate  bending  triangular  element  (called  discrete  Kirchhoff  triangle)  which 
allowed,  as  he  said,  a  significant  reduction  of  algebraic  operations  in  the 
evaluation  of  the  stiffness  matrices  and  bending  moments. 

In  1983,  Fricker  [Ref.  12],  based  on  Razzaque's  report,  described  a  new 
three-noded  triangular  element  for  plate  bending  with  an  extra  internal 
'bubble'  function  proposed  by  Irons  for  making  the  Razzaque's  element  more 
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flexible.  He  also  claimed  that  his  element  and  that  developed  by  Hansen, 
Bergan,  and  Syvertsen  were  the  two  most  accurate  triangular  thin  plate 
elements  currently  available.  Furthermore,  Lee  and  Zhang  [Ref.  13] 
developed  a  six-node  plate  bending  element  based  on  a  modified 
Hellinger-Reissener  principle  and  the  Reissner-Mindlin  plate  bending  theory 
to  obtain  more  accurate  and  reliable  solution. 

hi  this  thesis,  a  triangular  plate  bending  element,  which  was  developed 
by  Tocher[Ref.  IS],  is  used  to  analyze  the  laminated  composite  plates 
subjected  to  a  pure  bending  load  and  to  determine  the  deflections,  strains, 
and  stresses  of  the  plates. 

B.  TRIANGULAR  PLATE  BENDING  ELEMENT  FORMULATION 
The  plate,  in  general,  may  have  any  irregular  geometry  and  loading. 
The  derivations  are  based  cm  assumptions  of  the  small  deformation  theory, 
described  in  the  earlier  chapter.  A  triangular  plate  element  is  employed  for 
the  following  developement.  Each  nodal  displacement  of  the  element 
possesses  three  components:  a  displacement  in  the  z  direction,  w;  a  rotation 
about  the  x  axis,  wy  (derivative  of  w  with  respect  to  y);  and  a  rotation  about 
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the  y  axis,  wt  (derivative  of  w  with  respect  to  x),  as  shown  in  Figure  15.  The 

displacement  function,  w,  is  assumed  to  be 

*K*.y)  -  a \  +  a^x  +  a\y  +  a^x2  +  a$xy  +  a fjp  +  a-jx*  +  a%(x2y  +  xy2)  +  a 


-  M  a  . (33) 

where 

[n]  -  [1  x  y  x2  xy  f  x3  (xV+xy2)  y3  ]  . (34) 

and 

al  ' 
a2 

[«]-  ;  . os) 


i  i 

Here,  the  constants  a,,  a^,  .  ,  a,  must  be  determined  from  the  nodal 

conditions 

w(*.y)“q».  wy(x,y)-q2,  -ws(x,y)-q3  at  (x„y,)-(0,0) 

w(x,y)-q4,  wy(x,y)-q5,  -w^y)-^  at  (x2,y2)-(0,y2)  . (36) 

w(x,y)-q„  wy(x,y)-qs,  -wx(x,y)«q9  at  (x^-fx^ 

The  local  y  axis  is  taken  the  same  as  the  line  connecting  the  node  1  and  node 
2  with  the  origin  placed  at  node  1.  The  local  x  axis  is  taken  towards  node  3 
as  shown  in  Figure  16.  By  using  Eq.(33),  Eq.(36)  can  be  shown  in  the  matrix 


form  as 


Figure  15.  Nodal  degrees  of  freedom  of  a  triangular  plate 


in  bending 


Figure  16.  Local 


nodal  coordinate  system 


*9 


|(‘) 


(37) 


where  [n']  is  a  9x9  matrix  (see  Appendix).  By  using  Eqs.(33)  and  (37), 


Eq.(4)  can  be  expressed  as 

where 

0  0  0  2  0  0  6*  2y  0 
0000020  2x  6 y 
0  0  0  0  2  0  0  4 (x+y)  0 

and 


(38) 


(39) 


(40) 


Finally,  the  element  stiffness  matrix  in  the  local  (x,y)  coordinate  system  can 


be  written  as 


[H\TlQ)[H)dV 


(41) 


where  Vw  indicates  the  volume  of  the  element  and  the  matrix  [ Q ]  is  the 


material  property  matrix  for  the  generally  orthotropic  lamina.  By  substituting 


for  [H]  from  Eq.(40),  Eq.(41)  becomes 

[^e)]  -  (KT1 )  . (42> 

where  t  is  the  thickness  of  the  plate.  For  a  laminated  composite  plate 
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(43) 


-  }  Jf  (Cl*(4  -*J_, )  -  ID] 
and  the  components  within  the  bracket  of  Eq.(39)  won't  vary  with  the 
thickness.  The  element  stiffness  matrix  in  the  local  (x,y)  coordinate  can  be 
rewritten  as 


where 


[^e)]  •(["']■' 


lC]TlD][Cld4 


-1 


(44) 


[Cl- 


0002006*  2 y  0 
0000020  2*  6y 
0  0  0  0  2  0  0  4 (x+y)  0 


(45) 


The  area  integrals  appearing  inside  the  bracket  of  Eq.(44)  can  be  calculated 
in  terms  of  the  global  XY  coordinate  system  as  well  as  in  terms  of  local  xy 
system  chosen  in  Figure  16  from  the  following  relations: 


Sarea "  \x3  *2  (46) 

lxeadufym\Xly2  (47> 

gy^dxefy"  ^Jc3.y2  (y2  +^3)  (48) 

lj^.dcdy-±4y2  (49) 

§xy.  dxdy  -  +  2^3>  . (50) 

Jj\£  dxdy  -  ’^xyy2(y\  +yf)  . (51) 


Appendix  shows  the  detailed  oauculations  of  the  local  element  stiffness 


matrix.  Finally  the  element  stiffness  matrix  from  Eq.(44)  in  the  global 
coordinate  system  whose  XY  plane  is  assumed  to  be  same  as  the  local  xy 
plane)  can  be  obtained  from 

ra(e)-ip]7|'*<‘!>>i  . (52) 

where  the  transformation  matrix  [p]  is  given  by 

'10  000  000  O' 

0  lax  max  0  0  0  0  0  0 

0  loy  moy  0  0  0  0  0  0 

00  010  000  0 

[p]  “  0  0  0  0  lax  max  0  0  0  (53) 

^  0  0  0  0  loy  may  0  0  0 

00  000  010  0 

0  0  0  0  0  0  0  lax  max 

0  0  0  0  0  0  0  loy  may 

where  (lox,  m^)  and  (1^,  m^)  represent  the  direction  cosines  of  the  lines  ox 
and  oy  respectively. 

After  each  element  stiffness  matrix  is  calculated,  it  is  assembled  into 

the  global  stiffness  matrix.  The  resultant  matrix  equation  becomes 

WW  -  [F]  . (54) 

where  [F]  is  the  force  matrix. 

In  order  to  compute  stresses  and  strains,  post-processing  is 
accomplished.  After  obtaining  displacements,  local  strains  are  extracted 
from  the  global  displacements.  To  this  end,  the  local  coordinate,  matrix  [n']. 
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direction  consines,  and  the  transformation  function  for  each  element  are 
calculated.  Then,  the  local  strain  vector  is 

[Ve)  J  -  p/|V  j  lp][d\  . (55) 

where  matrix  [H']  is  the  same  as  matrix  [H]  except  that  x  is  one-third  of 
(Xj+Xj+Xj)  and  y  is  one-third  of  (y,+y2+y3)  and  z*+h  and  -h  in  order  to 
compute  strains  at  the  top  and  bottom  sides  of  the  plate.  In  addition,  [d]  is 
the  nodal  displacement  vector  of  an  element.  The  local  strain  is  transformed 
into  the  global  strain  by  using  the  following  equation: 

. (56) 

where 

lax  "*ax  lax*  max 

\T\m  nicx  lax  -lox*mox 

-2lox*max  2lax*max  lox~mox 

At  the  end,  stresses  are  computed  from  the  strains  using  Eq.  (6).  These  are 
the  global  stresses. 
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VL  NUMERICAL  RESULTS  AND  DISCUSSION 


A.  NUMERICAL  RESULTS 

In  order  to  calculate  the  strain  and  stress  for  a  laminated  composite 

plate  with  a  quarter  inch  hole  at  the  center  and  subjected  to  a  pure  bending 

moment,  a  program  named  MFEM.M  is  written  in  Matlab.  With  an  input  for 

material  properties,  plate  dimensions,  the  applied  bending  moment,  and  the 

failure  load  of  the  material,  this  program  calculates  the  deflection,  strain,  and 

failure  stress  for  each  element  Another  program  named  DMTRX.M  is  also 
written  in  Matlab  which  is  capable  to  calculate  the  stiffness  matrix  [Q]  for 
four  different  kinds  of  laminae  [0°],  [90°],  [0°/90°],2,  and  [(+45°/-45°)6]s. 

In  this  thesis,  a  laminated  plate  with  a  quarter  inch  hole  at  the  center  is 
modeled  with  214  nodes  and  371  elements,  which  generates  642  by  642 
stiffness  matrix  (642  dof)  as  shown  in  Figure  17.  After  calculations, 
normalized  stresses  of  stress  concentration  and  the  normalized  failure 
stresses  v.s  the  distance  away  from  the  circular  hole  edge  for  [0°]24  and 
[0°/90°]12  are  plotted  in  Figure  18  through  23.  As  a  result,  the  stress 
concentration  factor  calculated  from  the  F.E.  analyses  for  the  [0°]24  laminate 
is  1.47  and  that  for  the  [0°/90°],2  is  1.27.  Normalized  failure  stresses  for  the 
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[0°]24  laminate,  the  top  [0°]  layer  and  bottom  [90°]  layer  of  the  [0°/90°]l2 
laminate  are  all  close  to  one  at  the  edge  of  the  hole. 

B.  DISCUSSION 

Whitney  and  Nuismer  [Ref.  7]  said  that  when  the  stress  at  a 
characteristic  distance,  d,,,  away  from  the  discontinuity  would  be  equal  to  or 
greater  than  the  strength  of  the  unnotched  material,  failure  would  occur.  In 
this  current  work,  the  constant  value,  d^,,  for  the  [0°]^  and  [0°/90°]12  laminae 
is  zero  from  the  results.  Thus,  the  failure  criterion  of  Whitney  and  Nuismer 
assumption  is  not  suitable  in  this  kind  of  graphite  fiber/epoxy.  Since  it  is  of 
great  importance  for  design  purpose  to  be  able  to  predict  the  strength  of  a 
composite  under  loading  conditions  prevailing  in  service,  some  failure 
criteria  like  the  maximum  stress  criterion  and  the  Tsai-Hill  failure  criterion 
are  used  in  order  to  predict  the  failure  and  check  the  failure  criteria. 

1.  Maximum  Stress  Criterion 

This  criterion  tells  that  failure  will  occur  when  any  one  of  stress 
components  is  equal  to  or  greater  than  its  corresponding  allowable  strength. 
Since  the  stress  in  the  x-direction  is  the  dominant  stress  in  this  current  work, 
failure  will  occur  when 
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a  i  i.x\  . (57) 

where  XTt  is  the  ultimate  uniaxial  tensile  strength  in  the  fiber  direction  ( in 

this  case,  x-direction). 


2.  Tsai-Hill  Criterion 

According  to  the  Tsai-Hill  criterion,  failure  of  an  orthotropic  lamina 
will  occur  under  a  general  stress  state  when 


2  2  2 

gl  ql  a2  a2  ,  q12  >  , 

X 2  X2  Y2  S2 


(58) 


where  X,  Y,  and  S  are  the  longitudinal  tensile  failure  strength,  the  transverse 
tensile  failure  strength,  and  the  in-plane  shear  failure  strength,  respectively. 

Figures  7  through  14  show  that  for  this  particular  laminated  composite 
failure  occurred  almost  simultaneously  on  both  compressive  and  tensile 
sides.  From  figures  of  normalized  failure  stresses,  it  is  obvious  that  failure 
occurred  at  the  hole  edge  when  the  stresses  at  the  hole  edge  were  equal  to 
their  corresponding  allowable  strengths. 

Stress  concentration  factors  obtained  from  the  experiments  and 
numerical  analyses  are  nearly  the  same  for  both  [0°]24  and  [0°/90°]I2 
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laminated  plates.  This  result  confirm  that  failure  occurs  when  the  maximum 
stress  at  the  edge  of  a  hole  is  the  same  as  the  material  strength. 
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Figure  17.  x-y  plane  view  of  the  FEM  mesh  of  the  composite 
plate  with  a  circular  hole. 
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Normalized  Stress  01  S.C  F 


Figure  18.  Stress  concentration  factor  for  [01  laminate 
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Normalized  Stress  01  S.C.F 


1.3 


Figure  19.  Stress  concentration  factor  for  [0/90]  laminate  at  top  layer-[0] 
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Figure  20.  Stress  concentration  factor  for  [0/90]  laminate  at 
bottom  Iayer-[90] 
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Figure  21.  Stress  distributions  for  [o]  laminate 


Normalized  Failure  Stress 


Figure  22.  Stress  distributions  for  [0/90]  laminate  at  top  layer-[0] 
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Normalized  Failure  Stress 


Distance  Away  From  The  Hole  Edge  (y-fl)-(inches) 

Figure  23.  Stress  distributions  for  [0/90]  laminate  at  bottom  layer-[90] 


VI L  CONCLUSIONS 


In  this  thesis,  a  four-point  bending  test  is  used  to  investigate  the  failure 
strain  and  failure  stress  for  a  graphite- fiber/epoxy  composite  of  four  different 
laminated  plates.  From  experimental  results,  composite  elastic  modulii  are 
obtained  and  the  stress  concentration  factors  due  to  a  quarter  inch  hole 
located  at  the  center  of  the  plate  are  calculated.  Furthermore,  a  numerical 
analysis  using  the  finite  element  method  is  conducted  to  simulate  the  plate 
under  the  same  loading  condition.  By  using  both  experimental  and  numerical 
results,  the  characteristic  distance  d^,  which  determines  failure  based  on  the 
Whitney  and  Nuismer  theory,  is  examined.  Finally,  the  maximum  stress 
criterion  and  the  Tsai-Hill  failure  criterion  is  applied  to  check  the  failure  of 
composite.  These  failure  criteria  confirm  that  failure  occurs  when  the 
dominant  stress  component  of  the  laminated  composite  is  equal  to  its 
maximum  allowable  strength.  From  this  work,  it  is  also  concluded  that: 

♦  the  four-point  bending  test  is  a  good  way  to  investigate  material 
properties  although  it  is  a  little  more  complicated  to  conduct  than  the 
tensile  test. 
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♦  The  characteristic  distance  for  the  Whitney  and  Nuismer  theory  is 
zero  for  the  present  composite  material.  Both  the  maximum  stress 
criterion  and  Tsai-Hill  failure  criterion  are  proper  to  predict  the 
failure  for  this  material. 
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APPENDIX 


The  [n/]  matrix  in  Equation  (37)  is 
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The  detail  calculations  of  [C]T[DJ[C]  in  Equation  (44)  are  shown  as 
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where 


and 


la^Ditci 
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A4  *  (4£>]2  +  8D]g)x  +  (4D||  +8Dj5)y 
A5  -  (4D26  +  8 £>6$)x  +  (4D16  +  8 D66)y 
Ag  -  (4£>22  +  26 )*  +  (^12  +  ^ 26 

A7  -  (12jD|2  +  24Dig)x^  +  (12Dj j  +  2AD^)xy 

Ag  =  (4D22  +  16Z>26  +  ^6D66^  ^  ^D12  +  16£)16^ 

+(16£>26  +  32D66^  +  (4£>11  +  16Z)16  +  l6D66^y 


A9  *  (12Z>22  +  ^4^26^  +  (12/^12  +  24Z)26^ 

Thus,  the  area  integrals  inside  the  bracket  of  Eq.(44)  can  be  obtained 


$[c\Tmc\dA- 
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0000  0  0  0  0  0‘ 

0000  0  0  0  0  0 

0000  0  0  0  0  0 

0  0  0  i4j  A2  Ay  A 

0  0  0  i42  A$  A5  A%  A} 2  ^17 

0  0  0  A4  A$  Afi  Ag  A j3  -^28 

0  0  0  Ay  v4g  Ag  A-^q  ^24  ^19 
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where : 


Alm2Dux3y2 , 
A3*2D6fix3y2 , 


A2=2D,6x3y2 , 
A4=2D,2x3y2 , 
A4“2D2Jx3y2 , 


A7*,2D11x23y2 ,  A#m2D,6x23y2 , 

A9ax2D12x  3y2 ,  A10*3D,,x  jy2 , 

Au=(4D12+8DJ(x23y2)/6+(4Dn+8Di^[x3y2(y2+y3)/6]> 

A12“(4D26+8D66)(x23y2)/6+(4D16+8D66)[X3y2(y2+y3)/6]» 

A13*(4D22+8D25Xx23y2)/6+(4DI2+8D26)[X3y2(y2+y3)/6]» 

A14=(12D12+24DJ(x33y2)/12+(12D11+24D16)[x23y2(y2+2y3)/24], 

+(8D12+16D16+16D26+32D66)[x23y2(y2+2y3)/24] 

+(4Du+16D16+16DJ[x3y2(y22+y2y3+y23)/12]. 

Ai6s“2Di2x3y2(y2+y3),  Ai7“2D26x3y2(y2+y3), 

Ali-2D22X3y2(y2+y3),  Aig=3D12x  2(y2+2y3), 

A20=(12D22+24D26)[x23y2(y2+2y3)/24] 
+(12DI2+24D26)[x3y2(y22+y2y3+y23)/12], 

A^^D^  y2(y22+y2y3+y23). 
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